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We present an analysis of two-loop mixed QCD and electroweak corrections to 
the decay of the Z boson into light quarks. We find that the naive factorization 
of QCD and electroweak corrections does not describe correctly the two-loop 
effects. The nonfactorizable corrections shift the width of the Z boson by 
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approximately —0.55(3) MeV and increase the central value of the strong 
coupling constant determined at LEP by 0.001. 

With a total number of about four million hadronic Z decays collected at each of the four 
LEP experiments a final precision in the Z boson decay rate of about 2 MeV is expected. 
At the same time the ratio -Rhad — r(Z — > hadrons) /T(Z — > will be measured 

with an accuracy of about 0.1 per cent. These measurements will determine a s with a 
remaining uncertainty of about 5a s = ±3 x 10~ 3 . This puts not only severe constraints on 
the experimental analysis but also on the theoretical understanding of subtle higher order 
effects. QCD corrections have been calculated up to order a 3 |I| and leading as well as 
subleading electroweak corrections have been calculated up to two loops [0. 

Another important class of effects is provided by the interplay between electroweak in- 
teractions and QCD. The relation between Mw, Mz, Gf, and a is affected by the large top 
quark mass and again two and even three loop self-energies have been evaluated to arrive 
at accurate predictions. In addition there is of course the large class of "reducible" correc- 
tions which originate from one-loop electroweak diagrams, multiplied by the QCD correction 
factor (1 + a s /7t). The remaining two-loop effects are induced by nonfactorizable vertex cor- 
rections. With an order of magnitude of " wcak — fa 0.4 x 10 -3 and an unknown coefficient, 
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they are not a priori negligible for an analysis at the 0.1 per cent level. In particular it 
is evident that an ansatz assuming factorization between QCD and weak effects cannot be 
valid for the irreducible vertex diagrams. For the Z decay rate into bottom quarks this 
has been confirmed by the calculation of the leading contributions which are enhanced by 
a factor m^/rri^ and by logarithms ln(m^/m^) In both cases factorization is 

invalidated and nontrivial additional terms are obtained. 

These results were derived by expanding in the "small" parameter mf^/mf. In contrast, 



1 



no such expansion parameter is available for the Z decay rate into light u, d, s, or c quarks. 
All these quarks are effectively massless and both Z and W bosons appear as virtual particles 
in the relevant vertex diagrams which have to be evaluated at q 2 = m z . A similar problem 
arises also for the non-enhanced contributions to the bb vertex. However, the Z decay rate is 
dominated by the u, d, s, and c channels and the corresponding mixed terms will dominate 
in the total Z decay rate. For this reason we concentrate on light quark final state. 

The case of QCD combined with QED provides an illustrative example for mixed QCD- 
electroweak corrections. The factorization ansatz ^1 + Q^ff ) (l + ^r) i m P nes a mixed term 
Qoff However, the proper evaluation || leads to — fQqff which differs in magnitude 
and sign from the naive ansatz. 

With this motivation in mind we have evaluated the mixed corrections of order ° wcak — 

7T 7T 

for T u and Td which originate from irreducible vertex diagrams. Our approach is based on 
the observation that the relevant rates can be calculated for the decay of a virtual Z with 
mass squared equal q 2 , with q 2 alternatively far smaller or far larger than m z , by computing 
a large number of terms in the expansion in q 2 /m 2 z or m 2 z /q 2 . The results can then be 
extrapolated even to q 2 = m 2 z . 

The tree-level decay width of the Z boson into light quarks is given by 

K» ( z - „-) = + ( 3 ;f\ (i) 

9 q = Q q 

Cw 

(hq — S wQq) 
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where s\y and cw are sine and cosine of the weak mixing angle, I 3q and Q q denote the isospin 
and electric charge of the quark q, and iV^ = 3 is the number of colors. 

The Born level result (|l|) receives both QCD and electroweak corrections. The QCD 
corrections have been calculated up to three loops and yield a correction factor 

t qcd = p ( o) (2) 

First order electroweak corrections to the Z boson hadronic width can be divided up into 
three finite parts 

p (l loop EW) ( Z _^ g fl = T Z + T W + T ct (3) 

The first two contributions can be calculated from the imaginary parts of the on-shell self- 
energy diagrams shown in fig. ([!]). T z is given by the vertex correction (|l|a) with a Z boson 
exchange together with the Z boson contribution to the wave function renormalization of the 
quarks QT|t>) . T w is given by two analogous diagrams with the internal Z boson replaced by 
W, plus the diagram (JJc). In the linear 't Hooft-Feynman gauge, adopted in this paper, the 
sum of those three diagrams is ultraviolet divergent. It is made finite by including in T w the 
divergent part of the counterterm generated by the Z boson wave function renormalization; 
this is obtained by making the following replacements in the Zqq vertex: 
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(The calculation is done using dimensional regularization in 4 — 2e dimensions.) 

Finally, T ct is the tree-level decay width ([I]) multiplied by the remaining, finite Z boson 
wave function renormalization constant. 

The splitting of the decay rate into the three contributions is convenient for the descrip- 
tion of the QCD corrections. In particular, the effect of QCD corrections on T ct is just 
the multiplicative factor given in eq. (^). For the remaining two contributions to the decay 
Z — > qq we have 
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where q is the isospin partner of the quark q, s is the momentum squared of the external Z 
boson, and M is the mass of the virtual heavy boson inside the diagram. For a decay of an 
on-shell Z boson we have x = 1 in T z and x = M^/M^ w 1.292 in r . The functions A 2; 3 
can be expanded in a series in the strong coupling constant 
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where CV = 4/3 is the SU(3) color factor. The 0th order terms of these series have been 
calculated in |T0|,pJ] 
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The behavior of these functions at x — > and at x 

expansions for small and large external momentum. 

,(o) 



-> oo can be found from asymptotic 
It turns out that only the first few 
terms are needed to compute with good accuracy A\ u> for all values of the argument x. This 
observation is the basis of the present paper. Since virtual gluons do not radically change the 
analytical properties of the Z self-energy diagrams we can approximate the QCD correction 
functions A^ by the first few terms of their asymptotic expansion. In a large part of the 



computations we employ the program package MINCER |L2| written in FORM JT3] . 

The function A^ receives contributions from the diagrams shown in fig. (0) and fig. @. 
If the virtual particle is the Z boson the sum of these diagrams is finite. On the other hand, 
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in the case of the virtual W boson, the finite result is obtained only in the sum with the 
diagrams of fig. (||) and the counterterm discussed before eq. This is a consequence of 
the difference between the Z couplings to the quarks in diagrams |l|a and |I|b. 

We notice, however, that the diagrams in fig. (§) are proportional to {g~) 2 which can 
be rewritten as g q g^ + 2I^ q g~cw/sw- It is convenient to treat the two parts of this sum 
separately. The first part together with diagrams of fig. (fj) is finite and gives a contribution 
described by A2 . The second part can be combined with diagrams in fig. @ and the 
counterterm to give A3 . 

We can write down an expansion of A^ (x) around x = as a sum of three series 
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The first series is obtained from the finite part of the corrections to the heavy boson con- 
tribution to the quark wave function renormalization, fig. (|3|). The remaining two series 
describe the gluonic effects on the heavy boson vertex correction shown in fig. (0); it is 
convenient to isolate the terms containing C(3). The results are 
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In SY and S s we have not displayed terms O(x ) which are divergent but cancel in the sum. 

S s and Srf converge rapidly; in fact one can recognize the general formula for their 
coefficients and sum up both series exactly; the result (being a rather complicated function 
containing tetralogarithms) will be presented elsewhere. 

On the other hand, SY converges very slowly for x > 1; this is a consequence of a three- 
particle cut at s = M§ in the diagrams in fig. (fj). This cut corresponds to the decay channel 
Z — > W + du. If the internal heavy particle is a W boson we need the value of this series at 
x = 1.292. It is therefore necessary to compute the expansion of the function Ag on the 
other side of the three-particle cut, that is the asymptotic behavior at x — > oo. In analogy 
with eq. (§) we define 

A«(x) = T s {l/x) + lf(l/x) + T/(l/x)C(3) (12) 
In the sum T s + the divergences cancel and the first three terms give 

We do not present here the full formulas for the T-series for the following reason: beginning 
with the term 1/x 4 , the coefficients of l/x n are identical (up to the sign of the logs) to the 
coefficients of x n ~ 2 in the corresponding series S. Also in the terms 1/x 3 and x there is an 
equality in the sums T s +T± and S S + SY ■ This remarkable feature gives us confidence in the 
result. It should be stressed that the S and T series result from very different calculations. 

The equality of the coefficients guarantees that both expansions give equal results at 
x = 1. However, the slope of both approximations is not quite the same at x — 1; we can 
take the magnitude of the resulting cusp as an estimate of the numerical error in the final 
result. 

We also notice that values of A^ at x « 1, relevant for the mixed QCD/electroweak 
corrections, are well approximated by the asymptotic value Ag (x — > oo) = —3/8, which 
corresponds to the two loop mixed QCD/QED calculation ||. We use in the numerical 
analysis A^(l) « (1.292) « -0.37 ± 0.04. 

The calculation of A3 involves four three- loop diagrams shown in fig. (|4]). In the present 
analysis of the QCD corrections we neglect the influence of the real W emission [14"|-|16|1 . (We 



notice that a part of it, originating from W emission off quarks, has been included in the T 
series discussed above. In principle this partial treatment is not gauge invariant; the induced 
error is, however, negligible due to the phase space suppression of the W emission.) 

Taking into account only those cuts of diagrams in fig. @ which do not cut W lines we 
obtain 

a(D, n 1 L 89 2 116 3 53 4 

A\ '(x) = 14x H x 2 H x s H x 4 

3 w 72 V 60 525 1400 

+ 121275 X + 1009008 X + ° {X ] ) + 3x^ U 

29 3 1 1 

7 + — (1 + 21nx) + — (7-21nx) — (14) 

24 8x v ; 8x 2V ; 24x 3 v ' 
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The first part of this formula represents the finite part of the diagrams in fig. (f|); the 
remaining terms arise from the addition of the part of fig. (|]), as discussed above eq. 
We find A^(1.292) = -0.87(1). 

The mixed QCD/EW two-loop corrections are often naively estimated in the factorized 
form, in which one assumes that this effect is equal to the one-loop EW correction multiplied 
by the one-loop QCD factor a s /ir. Our study shows, however, that the assumption of 
factorization is misleading; for example the QCD correction to the function A2 has a relative 
minus sign. 

We can summarize our result as a difference between the full corrections and the factor- 
ization formula 
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The ambiguity of the finite part of the counterterm (eq. ^) cancels in this combination. 
Inserting into eq. ([TB] ) the values of functions A 
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and using a s = 0.12, a = 1/129, s^r 
of the nonfactorizable corrections is 



1.292) = 

0.223, M z = 91.19 GeV, we find that the net effect 
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The total change in the partial width YiZ 
down-type and 2 up-type quarks: 



hadrons) is obtained by summing over 2 
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-0.55(3) MeV 
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which translates into the change of the strong coupling constant determined at LEP 1 equal 
to 



AT(Z — > hadrons) 
T(Z — > hadrons) 
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(18) 



This shift is somewhat smaller but still of the same order of magnitude as the experimental 
accuracy and should to be taken into account in the final analysis of LEP 1 data. 
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FIGURES 




(d) 

FIG. 2. QCD corrections to the diagram 1(a) 




(g) 

FIG. 3. QCD corrections to the diagram 1(b) 
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FIG. 4. QCD corrections to the diagram 1(c) 
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